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Abstract. In this article, we construct loop soliton solutions and mixed soliton - 
loop soliton solution for the Degasperis-Procesi equation. To explore these solutions 
we adopt the procedure given by Matsuno [T|. By appropriately modifying the r- 
function given in the above paper we derive these solutions. We present the explicit 
form of one and two loop soliton solutions and mixed soliton - loop soliton solutions and 
investigate the interaction between (i) two loop soliton solutions in different parametric 
regimes and (ii) a loop soliton with a conventional soliton in detail. 
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1. Introduction 

In this paper, we investigate loop soliton solutions and their dynamics in the Degasperis- 
Procesi (DP) equation [2], 

Uf + 3k 3 u x — u xx t + Auu x = 3u x u xx + uu xxx , (1) 

where subscripts denote partial derivatives and k is a positive parameter. The 
integrability of the DP equation was proved by constructing a Lax pair, deriving an 
infinite sequence of conservation laws and the existence of a bi-Hamiltonian structure [3] . 
Eq.(PQ) arises in a hydrodynamical context jl]. Interestingly when k = 0, Eq. (pQ) admits 
peakon solutions which are of the form u = ce _ ' a_c *', where c is the velocity of the peakon. 
Consequently the N peakon solutions of the DP equation were constructed and the 
dynamics of these special solutions were also studied by Lundmark and Szmigielski [5] . 
The prolongation algebra and the Hamiltonian operator of this equation was reported 
in [6]. Vakhnenko and Parkes have studied various travelling wave solutions of the DP 
equation including hump-like, loop-like and coshoidal periodic-wave solutions j7]. Qiao 
has come-up with three new types of soliton solutions to this model, namely M-Shape 
peakons, dehisced solitons, double dehisced 1-peak solitons [8]. 

Here we are interested in exploring loop soliton solutions exhibited by Eq. (CQ). The 
motivation comes from the contemporary interest in studying loop soliton solutions in 
integrable nonlinear evolutionary equations P4T9]. To our knowledge periodic inverted 
loop solutions and one loop soliton solution were reported for the DP equation [7j. 
The explicit form of two loop soliton or higher order loop soliton solutions are not yet 
reported. Moreover the collision dynamics between the loop soliton solutions are also 
yet to be studied for this equation. 

To explore loop soliton solutions in Eq. ([[]) we follow the procedure given by 
Matsuno pQ. The author has derived iV-soliton solutions for the DP equation from the 
modified version of the Kaup equation since the r-function of the latter is already known. 
By carefully examining the r- function we observe that one has the freedom in changing 
the sign of the coefficients of the exponential functions so that even after the sign changes 
the resultant r-function satisfies the bilinear identities. By taking this advantage, we 
reshape the r-function appropriately and derive loop soliton solutions in a systematic 
manner. In the case of a one soliton solution, since there is only one exponential function 
present in the r-function we simply change the sign of the coefficient in front of this 
exponential function and obtain a one loop soliton solution. In the case of two soliton 
solutions since there are two exponential functions present in the r-function one has two 
choices in fixing the coefficients of these exponential functions. Either the coefficients of 
both the exponential functions are negative or one coefficient is negative and the other 
one is positive. 

For the first choice we obtain two inverted loop solitons and for the second choice we 
obtain a mixed loop soliton - smooth soliton type solution. We consider both solutions 
and study their interaction properties. To begin with, we allow two loop solitary waves 
(among the two, one is longer and the other is shorter) to travel in a particular direction. 
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As expected the taller loop wave travels faster than the smaller one and crosses the 
shorter loop wave in a finite time. When it approaches the smaller loop they start to 
interact. As a consequence the amplitude of the taller loop becomes shorter while the 
amplitude of the shorter loop becomes larger. This continues until the amplitudes of 
both the taller and shorter loop solitons become equal. 

In the second parametric regime, when the loop solitary waves cross each other 
the smaller and larger loop solitary waves overlap each other besides changing their 
amplitudes. In the third parametric regime we show that while the larger loop soliton 
overtakes the smaller one, the smaller loop solitary wave revolves circularly (in the 
clockwise direction) inside the larger loop solitary wave. 

In the second case (mixed loop soliton - smooth soliton) we allow the loop solitary 
wave to interact with the smooth solitary wave. Here we bring out a totally different 
kind of interaction. The smaller loop solitary wave travels along the surface of the 
smooth soliton. All these results are new. 

The plan of the paper is as follows. In Sec. 2, we recall the method of finding soliton 
solutions to this model. We present the method of constructing loop soliton solutions to 
this equation and derive the explicit expressions of one and two loop soliton solutions in 
Sec. 3. We investigate the collision dynamics, in the case of two loop solitons in detail. 
In Sec. 4, we derive mixed soliton - loop soliton solutions and study their interaction 
properties. We present our conclusions in Sec. 5. 

2. Method of finding soliton solutions [lj 

In this section, we briefly recall the method of constructing N- soliton solutions of (0Q). 
The DP equation can be written in a compact form, qi + q 2 u y = 0, first by defining a 
new variable q with 



and then introducing a reciprocal transformation, dy = qdx — qudt, dt = dt, in the 
resultant equation. Rewriting the relation (j2j) one can express the old variable (u) in 
terms of new variable (q) as 



where we have dropped the tilde above l V for simplicity. From (j3J) we can find u y . 
Substituting this derivative in the expression q t + q 2 u y = and differentiating the 
resultant equation with respect to y we obtain a third order ordinary differential equation 
in q. By successive differentiations this third order equation can be transformed to the 
first negative flow in the KK hierarchy (for more details one may refer [1]). The soliton 
solutions of this member can be identified from the soliton solutions of the modified 
version of the Kaup equation. From the known r-function of the Kaup equation one 
can go back and construct iV-soliton solutions for the DP equation. To do this first one 
should express q in terms of r-function of the Kaup equation. Doing so we get 



q 3 = u- u xx + ti 3 



(2) 



u = — 



q(lnq) ty + q 3 - k 3 , 



(3) 




(lnf) ty + K 2 , / = deU, 



(4) 
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where A = (a^) is a 2N x 2N matrix with elements 

a jk = (1 + S')8 jk + §^^-(1 - 
Pj -Qk 

isj-i = 6y = + c,-t - %o) + hi a,- 

P2J-1 = gj, <?2j-l = -Pj, P2j = Pj, 

with the soliton velocity defined by 



31 



j,k 

3 = 
hj ■ 



1,2, 



1,2,. 

= -Qs> 



N 
3 = 



1,2,. ..,N, 



3k 4 



1 - 



^k 2 



V 



K 2 kf 



3 



1,2 N. 



(5 a) 

(56) 
(5c) 



(6) 



In the above k/s are wave parameters and the amplitude parameters pj and q 3 - are found 
to be 



Pi 



1 + 



K/Jvj 




K 2 k"j . 



?i 



2 



1 - 




K2 7^2 



(7) 
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V 

- + 

K 



In a nutshell, to construct the solutions, one should first get q by substituting the explicit 
form of / (vide Eq. ( 15a}) ) in (j4|). By plugging this q in ([3]) one can get an implicit form 
of u(y,t). Finally, substituting the expression q in the mapping function, 

rv 1 1 

( )dy + d, (8) 

-oo q K 

where d is an integration constant, and integrating it one can get a relationship between 
x and y. We note that the function q(y, t) should not vanish for any y ant t. Otherwise 
a simple zero of q would yield a logarithmic singularity of x(y,t). Expressions ([3]) and 
(jHJ) constitute the solution for the DP equation in parametric form. Since q involves the 
variable y, after the integration, it is often difficult to express y in terms of x explicitly. 
In Ref. pQ, using the above procedure, the author has explicitly derived one and two 
soliton solutions of the DP Eq. (TjQ). 



3. Method of finding loop solitons 

To explore iV-loop soliton solutions of the DP equation we modify the bilinear solution 
(|5aj) as 

a jh = {l-(&)6 jk + ?^-(l-6 jk ), j,k = 1,2,. ..,N (9) 

with all other parameters pj, q~j and q k as given in Eqs. (j56|) -( l5~c|) . One can 
unambiguously prove the resultant r-function which comes out from the modified matrix 
elements also satisfies the bilinear identities. A simple sign change in the r-function leads 
to a new class of soliton solutions as we see below. 



3.1. One loop soliton solution 

Let us take AT = 1 in (j3J) with the matrix elements defined in Q. Expanding the 
determinant (jl]) we find the r-function modified to 

f = a 2 (l--e^+e 2 ^), (10) 
ai 
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where 

„ K 2 k 2 

3k 1 — 1 

£i = h(y+ ^^^ -yio), ai = ^ l _J k 2 , («*i>2). (11) 

One may observe that the sign in front of the second exponential in / is positive in the 
case of one smooth soliton solution [1] . However, this change of sign has other impacts 
on the solution. For example, in the regime nk\ < 2, q has a simple zero and u has a 
simple pole and the variable change y has a logarithmic singularity, as we see below. 

Substituting ({TO]) into (j3J) and carrying out the derivatives and simplifying the 
resultant equation, we find 

coshfi — 2ai + — 

= * ' ~ ai - (12) 

cosnfi 

Now plugging q and its derivatives in (j3]) and simplifying the resultant expressions to a 



compact form we arrive at 

3«» (ag-l)(a»-l) 

■u = — ^-p. (13) 

ax cosh^x — 2ai H 



ai 




Figure 1. Plots of (a) the inverted smooth solitary wave profile drawn u(y,t) versus 
y for the values k = 1.1, k\ — 2.1, (b) the one loop solitary wave drawn u(x, t) versus x 
with the same parametric values and (c) the inverse coordinate transformation between 
x versus y. 



Substituting q(y,t) in Eq.® and performing the integration we obtain the 
coordinate transformation between x and y in the following implicit form 

"1 + qji + (1 - «i)e ?1 



x 



V 

K 



In 



1 — ax + (1 + ax)e^ 



where 



(2ai-l)(ai + l) 



(14) 



(15) 



\ (2ai + l)(ai-l)' 

Since the variable y depends on £i one cannot express y in terms of x explicitly. 

Due to the coordinate transformation and the consequence of the change of sign in 
front of the exponential in the expression (fTOl) . the smooth solitary wave form morphed 
into a multivalued inverted loop like wave form with the condition given in Eq. fllip . 
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We present the nature of the obtained solutions in Fig. 1. To begin with, we draw the 
solution f fl3|) in the variable y where we get a smooth inverted solitary wave (Fig. la). 
While we depict the solution in terms of the original variable x we obtain only the loop 
solitary wave form (Fig. lb). To understand the geometrical connection between x and 
y we also draw a graph between these variables separately (Fig. lc). 



3.2. Two loop soliton solutions 

The r-function for the two loop soliton solution follows from Eq. (jl]) by restricting to 
N = 2. The associated determinant of / reads now, 

1 - aie fl 



/ 



Pl+91 

2pi 
P2+92 
Pi +Q2 
P2+g2 
P1+P2 



Pl+<?1 
2<?1 

1 - aie €l 

P2+<?2 
91+92 
P2+92 
P2+<?1 



£l±gl 

P2+qi 
Pi+gi 

P1+P2 

1 - d 2 e 6 
P2+g2 

2p 2 



pi+gi 

<Zl+<?2 
Pl+91 
Pl+<?2 
P2+g2 
292 

1 - a 2 e 6 



(16) 



When compared to two smooth soliton solutions the signs in front of the parameters 
ai and a 2 are different pQ. Substituting the expressions of p/s and q^s, j = 1,2, given 
in ([Zj), in the above determinant and expanding it we obtain the explicit form of the 
r-function as 

2u 2 



./ 



{0,10,2 

2 

Oi 



ai a 2 did 2 



,6+6 



e 

a 2 



26+6 



e 6+26 + e 2(6+6) 



(17) 



where 



v 



(fci-fc 2 ) 2 [K 2 (fc?-fcifc 2 + fcf) -3] 

(2Jb£ - A; 2 A; 2 + 2kj)n 2 - Q(k\ + fc 2 ) 



Substituting ffTTl) and (TT5|) in equation 
d(p, £), that is 

0102 



-3]' v ' 

we obtain a compact expression for 



q(y,t) 



9 

K— 
f 



f/(aia 2 ) 2 

The explicit forms of gi and p 2 are given by 
2 — c 2 — nk 2 



0i 



02 



5 



5 



2ai(l + nki) 

2 + Kfci 



2a 2 (l + Kk 2 
2 + kA; 2 



+ 



■e 6 + 



(2 - «Jfei)(2 - rejfe 2 ) 
4did 2 (l + Kki)(\ + k/c 2 ) 

(2 + ^0(2 + ^) 



,6+6 



,6+6 



(19) 



(20) 



(21) 



2ai(l — /cfci) 2d 2 (l — k/c 2 ) 4did 2 (l — ttfci)(l — k/c 2 ) 

We note here that the right hand side of Eq. (j3J) turns out to be a perfect square and the 
numerator in the resultant expression is factorized into the product of two functions, 
that is g — pip 2 with pi and p 2 being polynomials of e^ 1 and e^ 2 . 
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Substituting (|T9"|) in (E} and performing the differentiation we find the solution 
u(y, t) be of the form 

u = k 3 -, (22) 
9 

with 

9tfa 2 Aj gl 9<fc 2 A| > 9ft 2 fc 2 r gi+2£ 2 (oo\ 

ai (l-K 2 A;f) 2 a 2 (l - K^klf Q>\(\ — At 2 /c 2 ) 2 1 J 

a 2 (l - K 2 fc 2 2 ) 2 

We observe that the numerator of u is of the same form as for the two soliton solutions 
with the only difference in an over all sign change. The coordinate transformation 
between x and y is given by 

x (y,t) = V- + H-) + d, ( 24 ) 
where g\ and g<i are given in ( 120]) and ( 121]) . 




0.5 1 1.5 2 2.5 3 3.5 4 -3 -2.5 -2 -1.5 -1 -0.5 -8 -7.5 -7 -6.5 -6 -5.5 -5 -4.5 -4 



Figure 2. Two loop soliton solutions of (TTJ) drawn at different times, (a) t = — 10, 
(b) i = -0.8 and (c) f = 10 with k = 1.5, fci = 3.2, fc 2 = 3.8 and d = 0. 



5.5. Two Zoo;? soliton interactions 

Expressions (1221) and ( 124|) provide the complete description of the two loop soliton 
solution in the form of the parametric representation. 

It describes the two loop solitary wave troughs. Fig. (2) shows the interaction 
between two loop solitary wave troughs. Here we plot the solution u in terms of x 
with the parameter values k = 1.5, ki = 3.2, k 2 = 3.8 and d = 0. Both the solitary 
waves propagate towards the negative x direction. Initially at t = —10, the larger loop 
solitary wave is well separated from the smaller one (Fig. 2a). Since the larger loop 
soliton travels faster than the smaller one it starts to cross the smaller one in a finite 
time. Now the larger wave loses energy to the smaller wave. As a consequence the 
amplitude of the larger loop decreases and the amplitude of the smaller one increases. 
At t = —0.8, these two loop solitons have equal amplitude or energy (Fig. 2b). The 
two loop solitons do not superpose into a single wave as one can see from (Fig. 2b). 
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After the interaction (elastic interaction) these two solitons re-emerge and travel in their 
original direction by keeping their original amplitudes (Fig. 2c). 

Now we investigate the two loop soliton interaction in a different parametric regime, 
say for example k = 2.5, k\ = 3.4, k 2 = 4.8 and d = 0. To begin with (t = —5) the 
two loop solitons are well separated (Fig. 3a). As in the previous case both the larger 
and smaller loop waves move towards left. At t — — 1 the larger loop solitary wave 
approaches the smaller one (Fig. 3b). Fig. 3c shows the changes in their amplitudes 
when they start to interact with the amplitude of the taller one becoming shorter and 
the smaller one becoming larger. Interestingly at t — —0.25 the amplitudes of both 
the taller and smaller waves becomes equal and overlap each other (Fig. 3d). As time 
progresses the two loop solitons emerge out from each other in the original direction by 
changing their heights (Fig. 3e). At t = 5 the loops are completely separated and attain 
their original amplitudes. 

Another interesting loop soliton interaction can be seen by fixing the parametric 
values as k = 3.5, k\ = 10.4, k 2 = 4.2 and d = 0. We draw the interaction picture in 
Fig. 4. The two loop solitons are well separated from each other in the beginning, see 
Fig. 4a. When time moves on the larger loop meets the smaller loop (Fig. 4b) and 
they start to interact (Fig. 4c). In distinction to the previous two cases the smaller 
loop revolves inside the larger loop in the clockwise direction (Fig. 4c-4e) and finally 
emerges out from the larger loop. Later both the loops attain their original amplitudes. 
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Figure 4. The demonstration of clockwise rotation of a smaller loop soliton inside a 
bigger loop soliton. 



4. Mixed loop soliton - soliton solution of the DP equation 



In the previous section, while deriving the two loop soliton solutions we changed the 
sign of the coefficients of both exponential functions in the r-function to negative values, 
see Eq. ( fTBT) . In this section, we consider another possibility, we change the sign of one 
exponential function only 
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(25) 



For the opposite sign choice we would get inverted solutions. By performing the same 
procedure as in sub-section 3.1, we obtain the following expression for the functions /, 
q and u respectively, that is 
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(26) 
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where, 
9i : 



2 — nki 
2ai(l + k/ci 



-e ?1 + 



2 — k/c 2 
2a 2 (l + Kfc 2 ) 



-e>* - 
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4aia 2 (l + + nk 2 
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(29) 
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2ai(l — /tfci) 202(1 — ^2) 4aia2(l — — k^) 

ai(l - k 2 *; 2 ) 2 a 2 {\-^k\y ai(l — n 2 kf) 2 1 J 

9k 2 A| 2 £i+6 
a 2 (l - K 2 fc 2 2 ) 2 

Eq. (|2"5|) is a mixed soliton-loop soliton solution of the DP equation with the coordinate 
transformation given in Eq. (pM|) . 



4-1. Mixed soliton - loop soliton interaction 

In this sub-section we investigate the head on collision between a loop soliton with a 
smooth soliton. Eq. (128]) describes a smooth solitary wave form travelling towards right 
and a loop wave trough travelling towards left. These two waves propagate in opposite 
directions and collide with each other. To investigate the outcome we draw the solution 
( |28|) in Fig. (5) with the parameter values k = 0.91, k\ = 2.6 and k% = 0.91. In this 
case we observe that the anti-loop soliton creates a secondary wave crest. As time goes 
on the secondary pulse grows in amplitude while the primary wave loses its momentum. 
One can observe a double peaked smooth wave when the loop wave reaches the top of 
the smooth solitary wave. After a certain time the smaller loop starts sliding in the left 
direction. This interaction process continues until the two waves separate from each 
other. 
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Figure 5. The demonstration of mixed loop soliton versus smooth soliton interaction 
in ([!]). 
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5. Conclusion 

In this article, we have focused our attention on obtaining loop soliton solutions of 
the DP Eq. (Tfj). We have recovered these solutions directly from the r- function of 
the modified version of the Kaup equation. We have given parametric representations 
for both the pure loop solitons and mixed loop - smooth solitons and studied their 
wave dynamics. We have investigated the loop - loop soliton and mixed loop - soliton 
interactions in detail. In the case of the loop soliton interaction, when the amplitudes 
of the loop solitons are dissimilar, we observed that the smaller loop travels across the 
larger one in three different fashions before being shifted. We have shown the formation 
of double peaked soliton waves in the case of mixed loop - smooth soliton interaction. 
Currently, we are formulating N- loop soliton solutions for the DP equation. 
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